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Abstract 

We develop further the consequences of the irreducible-Boolean classifi- 
cation established in Ref. [9]; which have the advantage of allowing strong 
statistical calculations in disordered Boolean function models, such as the 
iVi^-Kauffman networks. We construct a ring-isomorphism {i ly . . . , i\} = 
V 2 [K] of the set of reducible i^-Boolean functions that are reducible in 
the Boolean arguments with indexes {ii, . . . , and the double power set 
V 2 [K] , of the first K natural numbers. This allows us, among other things, to 
calculate the number qk (A, (J) of i^-Boolean functions which are A-irreducible 
with weight to. qk (A, to) is a fundamental quantity in the study of the stabil- 
ity of A^i^-Kauffman networks against changes in their connections between 
their Boolean functions; as well as in the mean field study of their dynamics 
when Boolean irreducibility is taken into account. 
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1. Introduction 



Boolean functions play a seminal role in pure, and in applied mathematics. 
In pure mathematics: they constitute a field of study as themselves because 
of their rich structure. They appear in the theory of functional-graphs, as 
well as in the general theory of graphs, just to quote some related fields 1 . In 
applied mathematics: they support all the architecture of the modern digital 
computers as well as many fields of theoretical physics, theoretical biology, 
statistical mechanics, etc 2 ~ 5 . 

Our particular interest in them is mainly based in the work of Stuart Kauff- 
man which introduced the (now so-called) A^-Kauffman networks in order 
to try to understand how biological processes evolve; in an apparently spon- 
taneously way, from disorder to order 3,6 . iVi^-Kauffman networks are an 
example of disordered systems of many variables (Boolean functions) with 
deterministic rules of evolution. They are constructed in a random way, that 
trays to mimic the way in which nature constructs the basic bricks of the 
animated world; allowing a statistical treatment for them. One of the main 
problems with iVi^-Kauffman networks is to understand the way their dy- 
namics changes as the fundamental parameters of the model change i.e.: the 
number of Boolean functions N, its mean connectivity K, and its probabil- 
ity bias p that a Boolean function gives as output a "1" for a preassigned 
input 7-11 . This, however has only been done in an exact statistical way 
(including asymptotic expansions), only for special values of N, K, and p: 
The case of the so-called random map model where K = N, and p — 1/2 (so 
the Boolean functions are extracted with equiprobability) 7 . And the case 
K = l, also with p = 1/2 8 . 

In Ref . [9] , a new classification of Boolean functions, in terms of their real 
connectivity, was established as a necessity to understand new aspects of their 
dynamics as a statistical system. This led to the concept of irreducible- degree 
of a Boolean function, and allowed to understand new asymptotic properties 
of ATif-Kauffman networks 9-11 . We stress that this classification should be 
not confused with the one developed by Kauffman, which divides Boolean 
functions into canalizing and no-canalizing 3 ' 12 . While canalizing classifica- 
tion is a first step to understand the dynamics of AiT-Kauffman networks; it 
gives less information about a Boolean function, than its irreducible-degree. 
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As examples: i) The injectivity properties of the genotype-phenotype func- 
tion \I/ which is responsible of the dynamical diversity in iVX-Kauffman net- 
works, are obtained by an expansion in terms of the irreducible-degree 9,11 . ii) 
The calculation for the probability that an iVA-Kauffman network remains 
invariant against a change in one of their ^-connection functions depends 
on a series that involves term by term the irreducible-degree of the Boolean 
functions 9 . See also Eq. (Al). Hi) The dynamics of AA-Kauffman networks 
study through the mean field approach is based in their mean connectiv- 
ity, which is regulated by the irreducible-degree, and not by the canalizing 
classification. See Appendix A for a detailed discussion about canalization. 

While our specific research is based in Kauffman's networks, the need of 
new mathematical tools to deal with them has shove us to study the rela- 
tion of Boolean functions with rings. Thus, in this work we construct rings 
for the sets of Boolean functions which are reducible in some of their argu- 
ments and use this powerful isomorphism to calculate the number of Boolean 
functions with irreducible-degree A, and weight u, denoted by qk (A, a;). We 
also established an analytic expression for the irreducible-degree A of a given 
Boolean function. 

The content of this article is as follows: In Sec. 2, we introduce the formal 
apparatus, defining the concept of irreducible-degree of a Boolean function 
originally proposed in Ref. [9]. We establish their classification in terms of 
it, and some previous developed counting formulas are settled. It is also 
introduced the concept of weight, and some useful decompositions of the 
space of Boolean functions are shown. In Sec. 3, we recast the terminology 
of Boolean functions into the language of set theory by associating to each 
Boolean function on K arguments a set B which is an element of V 2 [K] , the 
double power set of the first K natural numbers. With this framework we 
further study the mathematical structure of the Boolean functions. Sec. 4, 
contains our main results; there we find: i) An analytical expression for 
the irreducible-degree A (bx) of a Boolean function b K on K arguments. 
ii) A subring of the set of reducible functions on A arguments, of the ring 
(V 2 [A], A, n). Hi) Last but not least, we find an analytic expression for 
g(X,u), the number of Boolean functions bx with irreducible-degree A and 
weight oo. In Sec. 5 we discuss our results. In Appendix A, the mathematical 
definition of a canalizing function is introduced, and the difference with their 
irreducible-degree is explained in detail. In Appendix B we make a summary 



3 



of the treatment developed in Ref. [10], where a study of the phase transition 
curve of iVif-Kauffman networks is done. 



2. Boolean Functions and their Irreducible Classification 



Along this work: We use the convention that the set of natural numbers 
starts with the number 1, i.e. N = {1, 2, 3, . . . }. And we denote by [K] = 
{1,2, ... ,K} the set of the first K natural numbers. For all S,S G Z 2 , 
additions S + S are intended to be modulo 2. When clarity will be required 

we emphasize this by S + 5 J . 

Definition 1: 

i) A K-Boolean function is a map 

b K : Zf ->• Z 2 , (1) 

with the assignment 

S = (S 1 ,...,S K )^a. (2) 

m) Its negation, denoted by -\>k is given by 

-nb K (S) = b K (S) + l. VSGZf (3) 

iV.5. There is a total order among the inputs S G Zf of 6^; which we are 
going to follow along this work. It is given by the bijection Zf- ^ — > [2 K ] , 
given by 

K 

s (S) = 1 + S i T ~ l 1 < s (S) < 2 K . (4) 

i=i 

Definition 2: 

A i^-Boolean function (1) is completely determined by its truth table 
tB(b K ), given by 

*B (b K ) = [0-1,0-2, ... ,<t 2 k\ , (5) 
where, a s G Z 2 , is the s-th image of (2), under the total order (4). 
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There are 2 2 ^ If -truth tables 05(6^), (5) corresponding to the possible 
Boolean functions (1). Thus 23 defines a bijection and we make the 

Definition 3: 

The set of all K -Boolean functions is given by 

~ K = {b K : Zf — ► Z 2 } A Zf . 

It can also be given a total order to the elements of E K by the bijection 



Zf A. 



2 2 " 



, defined by 



// = 1 + ]T 2 s - V s 1 < /i < 2 2 " , (6) 

s=l 

which is usually called Wolfram's classification 5 ; and is the one used in 
Refs. [9,10,11]. 

An important function associated to the elements of S# is the weight func- 
tion oj :E K — > [2*"] U {0} defined by 



(b K ) = (7) 



s=l 



with their values denoted by u. The weight function appears, mainly, when 
stochastic extraction of the functions bK is involved, as is the case with NK- 
Kauffman networks. There, each a s G Z 2 of the truth table (5) is extracted 
with a bias probability p (0 < p < 1) that a s — 1, so bx has a probability 

n(b K ) = n(uj)=p» (i-p) 2K -" (8) 

of being extracted 9 ' 10 . Due to this important application the weight function 
is going to be considered further in the work. 

We clarify our notation by the example of the truth table (5), for the 
case K = 2; presented in Table 1. In the first column, the images a s of 
the Boolean functions are arranged according to the total order of their two 
inputs (Si, S 2 ) given by (4). The next columns enumerate them according to 
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their Wolfram's number (6), which is indicated by numbers in boldface. At 
the bottom of the Table, 5 stands for the logical meaning of each 2-Boolean 
function. For \x = 11,13, $ — k — Si (i — 1,2) stands for the identity 
2-Boolean function in the i-th argument, while for /j, — 4, 6, $ — = ->Si 
represent their negations (3). The parameter A, of each Boolean function, to 
be defined shortly is its degree of irreducibility; while u, its weight; given by 
(7). 
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Table 1. The 23 (62) truth tables of the sixteen 2-Boolean functions. 

Not all the i^-Boolean functions depend completely on their K arguments. 
Extreme examples, for K = 2, from Table 1: are rules 1 and 16 (contra- 
diction and tautology, respectively), which do not depend on either S 1 or 
S*2. Rules 4, 6, 11 and 13 only depend on one of the arguments; while the 
remaining 10 depend on both. It has been shown that these facts are re- 
sponsible of important transition behaviors in A^-Kauffman networks, with 
repercussions in their applications to biology 9-11 . Then, a classification in 
terms of what is going to be called the irreducible- degree of bk] was proposed 
in Ref. [9]. 

Let us make the following 
Definitions 4: 

i) V % e [K], b K is {i}-irreducible <=> 3 S = {S u ...,S K )e Zf : 

bx (Si, . . . , Si, . . . , Sk) = 1 + bx (Si, . . . , Si + 1, . . . , Sk) ■ 
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On the contrary, if V S G Zf 

bx (Si, . . . ,Si, . . . , Sk) = bx (Si, . . . , Si + 1, . . . , Sk) ] 
bx is {i}- reducible. 

ii) V {ii, . . . , i\} C [X], 6^ is {ii, . . . , irreducible if it is {i a }- irreducible 
for a = 1, . . . , A. Similarly, bx is {i±, . . . , i\}- reducible. 

Hi) bx is irreducible of degree A (A = 0, 1, . . . , K); if it is irreducible over 
A indexes . . . C [K] and reducible on the remaining K — A 
indexes [K] \ {ii, . . . , 

iv) For \ = K, bx is called to be totally-irreducible. 

v) V i G [X] 

£Hx {i} = {6x £ | 6x is {i} — reducible} 

vi) V {ii,...,i x }<Z[K] 

A 

9Mii,...,i A }= n*K*:{i a } (9) 

a=l 

N.B. Vlx {h, ■ ■ ■ carries no information about whether, their ele- 
ments, are (or not) reducible on the remaining . . . , i\} indexes. 

vii) V bx G B.k, A (bx) represents the function 

\:~x^[K]U {0} (10) 
that gives the irreducible- degree of b K , with A standing for their values. 

From these Definitions we have the following 

Decompositions: 
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i) E K is disjointed decomposed in terms of the irreducible-degree (6x) of 
their elements by 



K 



A=0 



where 



1 K (\) = {b K e~K\\(b K ) = \}. 



(11a) 



(116) 



ii) From (9), %k (A) may be decomposed in its turn, as a disjoint union 
over all the indexes {ii, . . . , i\} C [K] with cardinality A, in the follow- 
ing way: 

a) The i^-Boolean functions . . . irreducible are given, from Defi- 
nitions 4 (H) and (vi), by 

a / A 

f| (e k \ m K {t a }) = e k \1\Jx k {i a } 



a=l 



where Morgan's Law was used. 



b) The fT-Boolean functions which are reducible in j G [K\ \ {ii, . . . ,i\} 
are given by f]j <R K {j}. 

Then 



(A) = [J 

{h,...,i x }C[K} 



n 



n 



n 

je[K]\{h,...,i x } 



(12) 



Hi) Ek is disjointed decomposed in terms of the weight u (bx) of their 
elements by 

2 K 

|J ¥k (w) , (13a) 



w=0 



where 



(u) = {b K eE K \u (b K ) = uj} , 



with cardinalities 



(136) 
(14) 
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A recursive formula for the cardinalities (5k (A) = #1^ (A) was found in 
Ref. [9] obtaining 



Pk{X)= {\) ® A ' (15a) 
where <3\ = (5\ (A). Taking cardinalities in decomposition (11), it follows, 

K 



A=0 

which can be inverted to obtain 13 : 

A 



0a = £(-1)W A J 2 2 ™. 

m=0 ^ ' 



(156) 



All the coefficients (3 K (A), but (3 K (0) = 2, grow with K. Note that % K (0) = 
{-it, t} consists only in the contradiction, and tautology functions; and their 
truth tables (5) are given by 

9M-.r) = [0,0,...,0], 

and 

^ K {r) = [l,l,. ..,!]. 

V v ' 

2 K 

On the other hand, from (15b) an asymptotic expression, for the number 
of totally- irreducible functions (3k{K) = <5k, for K ^> 1, is obtained with 
respect to the total number of if-Boolean functions 2 2 , giving 

g«i-o(^). <i«) 

So, with respect to the normalized counting measure, almost any fT-Boolean 
function is totally-irreducible. 

Irreducibility in Boolean functions shows us, that the real connectivity 
of a i^-Boolean function bx is not K; but A (bx)- So f3x(ty, defined by 
(15), gives a real gauge for it. When stochastic extraction of bx through (8) 
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is involved, the function uj (bx) given by (7) norms the average amount of 
Boolean functions involved in such processes. So, an important quantity to 
be considered when both reducibility, and stochastic extraction of bx must 
be taken into account; is the joint probability distribution in terms of the 
irreducible-degree, and the weight, given by qk (A, a;) /2 2 , where 

q k {\oj) =#[Xx(A)n<£*: {u)\. 

Calculation of qk (A,cu) is quite involved if strict use of Decompositions (11) 
and (13), with their cardinalities (14) and (15) are only used. In the follow- 
ing sections, we will construct a ring structure support to deal easily with 
quantities like this one. 

3. fT-Boolean Functions and their representation in V 2 [K\. 

It is possible to obtain a better understanding of Boolean irreducibility 
by recasting the description from Definitions 4, into set and ring theoreti- 
cal languages; which will increase considerably our calculation combinatorial 
counting power of important quantities. More concrete, (2^, +, •) consti- 
tutes a Boolean ring, and we are going to find subrings related to sets of 
reducible functions, to be defined below. Let us see this in detail: 

In general, for any set Q, its power set V Q = {T> \ T> C Q} is a Boolean 
ring with the set operations symmetrical difference A (addition), intersection 
fl (product), and; with and Q constituting the identical elements under 
addition and product, respectively 14 . There is a ring-isomorphism <£> into 
the set of Boolean functions 

2 n = {X:Q^Z 2 } A?S], (17) 

by the assignment of a characteristic set T>x of the boolean function X; which 
is defined by 

V x = $ (X) = {d E n | X (d) = 1} C fi; 

and has as a unique inverse association the characteristic function Xv of the 
set V, defined by 

r 1 ifdeV 

X v (d) = I vvevn . 

(O ifd<£V 
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So, VD^ePfi, the following ring-isomorphic properties are satisfied 
i) Xt>at — Xp + 

H) y.T>C\T = X© ' Xjr 

Hi) X$ = 0, Xq = 1. 
See details in Ref. [14]. 

Using the ring-isomorphism (17), for the inputs S G Zf of (2), we have 
z k ^ 2 [K] = | X . ^] _^ z ^ p ^] 

with the characteristic set Ax of X given by 

A x = r (X) = {i G [X] | X (i) = $ = 1} C [X] , (18) 

where Si is given by (2). The bijection (4) establishes a total order s 
(s — 1, . . . , 2 K ) among the elements of V [if]. So, we may label them by 

As = c [X], 

where: i Q G [if], a G [/], and / G [if]. Tautologically, 

K 

8 = s(A) = l+Y i T*- 1 = 1 + 2 <_1 Xa (0 • (19) 

Same consideration may now be applied to S^, with the total order (19); 
by using once again the ring-isomorphism (17). We obtain; 

E K = 2^1 = {X : V [if] ->• Z 2 } A [if] = [if] , (20) 

with the characteristic set Ex of X given by 

B x = V(3t) = {A,eV[K] \ X(A,) = <t 8 = 1}CV[K], (21) 

with <r s given by the truth table (5) of the if-Boolean function X. Now, 
using bijection (6), we may use the label \i (// = 1, . . . , 2 2 ) for the elements 
of V 2 [if], to obtain: 

B, = {A S1 ,...,A S JCV[K], 
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where: sp G [2 K ], (3 G [m], and m G [2 K ] . Tautologically 

2 K 

fi = fi{B)=1+ 2- 1 = 1 + ^2 S - 1 X B (A s ). 

s€{s\A s eB} s=l 

The following associations are going to be done in the future: 

B = * (6k) g [X] , 

and through (18) and (19) 

A = T(S) EV[K]. (22) 
To be noted also from (7) and (21) that 

u>(b K ) = #*(b K ) = #B. (23) 



4. fR K {h,---,ix} is a Subring of (V 2 [K] , A, n). 



For 91^ {ii, . . . , given by (9), we are going to show that there is a 
ring-isomorphism D\ K . . . , i\} = V 2 ([K] \ . . . , that allows us to 
easy count weight functions by means of (23). To do so, first we will recast 
Definitions 4 for the elements of E K in the language of the elements of V 2 [K\. 
Let us begin with the following: 
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Lemmas 



i.a) V^eP[if],andVi G [if] A ^ A A {0, and ( A A {i} ) A {i} = A. 

Proof: Follows from the fact that {i} ^ and the nilpotent property 
of the symmetrical difference A. 

■ 

i.b) V A G V [K], and V i, j G [K], such that i^j 
AA{i}A{j} = AA{i,j}. 

Proof: Follows from the fact that {i} fl {j} = 0. 

■ 

ii) Let bx '■ 1% — > Z 2 , and B = * (&#-) G "P 2 [if] be its associated set. 
Then: b K is {ij-reducible <^ V A G 7? [if], X B (A) = X B (A A {i}) . 

Proof: From (22), and Lema l.a; 

AA{t} = T(S 1 ,...,S l ,...,S K ) Ar(J_41,0,..,0) 

i 

= r (s*!, . . . , Si + i, . . . , Sk) ■ 

Now, from Definition 4-i the Lemma follows. 
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Corollary I 

Vb K e~ K \{b K ) = \^b K ). 

Proof: Let B = (-.&*), from (21) X B (A) = X B (A) + l, V A G V [K]. 
So, from Lemma (ii); bx is {i}-reducible •<=>- —16^ is {i}-reducible. 

■ 

Lemma 

Hi) With the total order of V\K\ given by (4), or equivalently (19), 
VAeV[K] 

s {A) = s{AA {t}) + T- 1 (2 X A (i) - 1) . (24) 
Proof: V A G V [K] , the inverse image of (22) gives 

r- 1 (A) = (s 1 ,...,s i ,...,s K ), 

and 

r- 1 (AA{i}) = (S 1 ,...,S i + l,...,S K ). 
Then from (19) 

K 

s (AA{i}) = l + J2 s i + + ^2 

= 1 + 5] S j V- 1 -2 i - 1 (S i - [Si + 1] 2 ) 
3=1 

= s(A)-2 i ~ 1 (2X A (i)-l). 

■ 

w) V£ = * (6 X ) G [K], with the total order (19), and with the truth 
table 23 (b K ) given by (5): 

b K is {2} - reducible \M G T 3 [K\ a s{A ) = cr s ( A A{i}) • (25) 
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Proof: From Lemma (ii), and the ring-isomorphic association (20): 
(J 8 (A) = X B (-4), and (T s (AA{i}) = %b (A A {«}). 



Theorem I: The operational form of the irreducible-degree (10) of 
b K G E K is given by 

K 

\(b K ) = J2 ©(-MM), (26a) 
i=i 

where, for any a G R, the step function 6 is given by 

1 if a > 

6 (a) = { , (266) 

if a < 

2 K-i (fe-l)2 l +2 i - 1 

T K ip K ; i) = ^2 [<?s + ^+2*-!)] 2 • (26c) 

6=1 s=(6-l)2'+l 

Proof: From (24), and (25) follows that if bx is {i}- irreducible, there 
exists an A G V [K] such that <7 s (a) 7^ °"s(^A{i})- Arranging sets, from 
(24) through the total order (19), so that, A ~< AA {i}, whenever 
s(A) < s{AA{i}) it follows that s(A) = s{AA{i}) - 2 i_1 : So, 
the double sum checks over the total number of 2 K ~ l pairs of indexes 
[s (A) , s (A A {«})], if it happens that cr s („4) 7^ o" s (.4A{i})- If that is 
the case at least for one of such pairs; then Tk (Pk', i) > 0, and so 
(J-r (pK'i i) ) — 1; implying that is {i}-irreducible. 



and 



Lemma 

u) V i G [if] there is a ring-isomorphism 

P 2 ([if]\{0) -^W* {iJC? 2 [K] 

given by 
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where 

B, t = {AA{i} eV[K] | A e B} , (27) 
and the unique inverse association 

<pr 1 (B) = {AeB\i(£A}. (28) 

Proof: From Lemma (ii); Vi (B) G m K {i} for all B G V 2 ([K] \ {i}). 
So, due to the unique inverse association (28) ifii is a bijection. Now, 

from (27) follows directly that V B,B eV 2 {[K\ \ {%)) A#) = 

BiABi, and (b H Bj = Bi n B { . Since, by (27) and Lemma (i.a) 
B fl = 0, it follows that £> U #j = 23 A <Bj. Then the ring operations 
A and fl are easily handled to show that ipi (b A Bj = ipi (B) A (fi(B), 

and ^ (B n B) =<Pi(B) n ^(B),\/ B,BeV 2 ([K]\{i}). 

■ 

Corollary II 



mK {i} = #V 2 ([K]\{i}) = 2 2 



Theorem II: V {ix, . . . , i\} G V [K] there is a ring-isomorphism 
: V 2 ([K] \ {h, . . . , * A }) — ► 9^ {ij, . . . , i x } C [X] 

given by 

/ \ 

j J ^ia,i( 



<p h _ h (B)=Bu[\jB la \u 

Ma} 



u . . . 



y{«a 1 ,«Q, 2 } y 



u 



I I Bi a ^..^i a . 



u • • • u B. 



(29) 
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where [i ai , . . . , i aj ] C («i < ••• < a,-) runs over all the 

subsets of {ii, . . . , i\} in cardinal order j — 1, . . . , A, and, 

B iai _ ia . = {AA{i ai ,...,i aj }eV[K] \ Ae B} . (30) 

With the inverse given by 

<P7*.., ix (B) = {AeB\i a tA Va = l,...,A}. 

Proof: From (30) and Lemmas (i), V B G V 2 {[K] \ i x }) 

BHBi i =Bi i n B in iR = 0, 

V {i Ql , . . . , i aj } 7^ {i^, . . . , i/3 fc }- So the union is disjoint. By induction 
from Lemma (v), aided by Lemma (i.b), and the property that 

Be9tA-{ii,...,iA} 

[ V ^4 G i3 -<=>- A A C <E B; VCC{i!,...,i A }]; 
which follows from Lemmas (i.b) and (m,); the Theorem is proved. 



Corollary III 



{ii, . . . , < A } = #P 2 ([X] \ {ix, . . . , *a}) = 2 2 * 
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Lemma 



vi) There is a bijection 



X K {i 1 ,...,i x }ny K (uj)^{BeV 2 ([K}\{i 1 ,...,i x }) = . 

Proof: From (30); = #B iai ,...,i aj V {/,, /,, } C {h,...,i x }- 

Taking cardinalities in the disjoint union (29) we have 



#^,..., iA (B) = #BJ2 (fi = * B 2A • 

3=0 V<7/ 

So { contracts cardinalities, and from (23), the weights by a factor 2~ A . 
From (13b) the Lemma follows. 

■ 

Corollary IV 



CO 



)K-\ 



where, V a e 



f 1 ifa = 

o(a) = 

[ ifa^O 

is Kronecker's delta, and [a\ e Z the /?oor function, defined as the greatest 
integer [a\ ; such that [aj < a. 

Proof: For any set B it happens that |BgN U {0}. From Lemma 
(vi): 

VKk {ii,...,ix}ny K (w) = 
whenever w/2 J pU {0}. Otherwise V6e^ {Ji, . . . , i A } n ^ (w) 
3 6 C V([K] \{h,...,ix}), such that #B = oo/2 x with B = 
<Pii,...,ix (&)■ 
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Theorem III: g K (A, u) — # [ %k (A) fl tyx (w) ] is analytically given 
by: 



A 



2 m 

x I M ) 

, L 2 ^-™ J 



x 



2^- 



(31) 



where 0° = 1. 



Proof: Taking the intersection of (12) and (13b), and using the idem- 
potent property of set's intersection (V fl V = V for any set V) we 
obtain 



{i u ...,i x }C[K] 



Zk\[{JXk {la} 



\a=l 



n 



u 

{il,...,*A}C[X] 



j6[X]\{ii,-,ix} 



^ H \ |J m K {i a } n (w) 



va=l 



n 



n 



f| «K*{j}n<p*(w) 

J6[K]\{<1.-.*A} 



u 

{il,...,i A }C[jq 



f| ^0'}n^H 

je[K]\{h,...,i x } 



U^HniRHU n f| JH^O'} 

je[x]\{n,...,i A } 
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u {[!u\ 2 :. »}■ 



(32a) 



{»i,...,i A }C[tf] 



From Theorem II 



je[^]\{n,...,i A } 



from Lemma (vi) 

Vk H n .w,- : , A = (bgp 2 ({<!, . . . , ^ a}) I = ^ 

and from Corollary IV; 

2 A 



#[iL 



U,...,lA 



2^- 



2 K- 



(326) 



Since, [2]^ ^ is a union of, non necessarily, disjoint sets; we have 



14 



#[2] 



— >»A 



E(- 1 ) n_1 



X 



n=l 



3=1 



(32c) 



Now 



j6([Jsr]\{u,...,n})u{ja 1 ,...,i«„} 

fl 0'} = ^ (K • • • , n} \ |~V , • • • , iaj) , 

je[/il\({il,-,»A}\{iai,-,»a„}) 

where Theorem II has been used; and it is worthwhile to bear in mind that 
{i ai , ■ ■ ■ ,j a J ^ Oi, • • • From Lemma (vi) and Corollary IV we obtain 

2 A-n 



# 



6£;:::!L. n (<") 



|_ 2 K ^A+nJ 
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2K-\+n 



Going to (32a), (32b), (32c), and taking into account that the cardinal count- 
ing result is independent of . . . C [K]; an overall factor (^) is ob- 
tained in the union (32a). 



Checks: The following formulas, for the number of A-irreducible func- 
tions with weight u, come as a result of consistency of (31) with (14), 
and (15): 



^2 Qk(\w)={ J and g K (A, u) = (3 K (A) . 

A=0 \ U / w=0 

See Appendix B of Ref. [10] for details on calculations. 
Special Values and Properties: 

1) 

q k (0,uj) = 5(uj) + 5(uj-2 k ), 

and 

Qk(\0) = Qk (\,2 k )=5(X). 

That is; the only completely reducible i^-Boolean functions are the 
contradiction ->t and tautology r functions, and also are the only two 
that have the most extreme values of cu: and 2 K respectively. 

2) 

g K (l,co) = 2K 5 (to -2 K - 1 ) . 

Which corresponds to the K-identities Li = (Si, . . . , Sk) = Si as- 
sociated to each one of the arguments Si G Z 2 , i = 1, . . . , K; and the 
corresponding K negations -itj = -\>k (Si, . . . , Sk) = Si + 1. 

3) qk (A, o>) is a symmetrical function of u at the value 2 K ~ 1 i.e.: 

q k (\,u) = q k (\,2 k -u). (33) 
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Proof: From the definition (7) of oj(b K ), uj{-b K ) = 2 K — u(b K ). From 
Corollary I, A(-ife^) = Xi^bx), then: V bx £ ^k- 

b K e % K (A) n <$k M e (A) n (2* - u) 



4) For w = 2n- l,n G N, 



2 



^(A,2n-1)= (^ 2n _ 1 J 5(^-A). (34) 
Proof: Follows directly from (31). 

■ 

5) For w = 2n, n E N, and if > 1, 

^(if,2n)« [l-2l(if,n)] 

where the function 21 (if, n) goes to zero faster than o (^) for n ~ 
(9(1) as if grows; and faster than o ( 22 k-i ) for, n ~ 2 A: " 2 (in the 

region of the maximum). So, for if >> 1, qk (if, 2n) obeys a Gaussian 
probability distribution with the same moments as for the odd case 
(34). 

Proof: From (31) with A = if the leading term is Of the re- 

maining if terms the next one in size is ( 2 ) . Using Stirling's ap- 
proximation for the factorials 



2*"^ ^ (2 K 
n ) ' \2n 



o(£r) forn~C(l) 

o{^k=r) forn~2 K " 2 
each one of the remaining if — 1 terms giving smaller contributions. 

I 

Corollary V 
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Let 



3 K = {b K EE K \co (b K ) =2n-l,neN}, 



then 



3k £ Ik (K) . 



Proof: Follows directly from Theorem III. 

That is, any i^-Boolean function bx with an odd weight is totally-irreducible. 
N.B. The converse is not true, for example, in the case K — 2 of Table 1, the 
2-Boolean functions 7 = and 10 = both have w = 2 and are totally- 
irreducible; while the other eight totally-irreducible 2-Boolean functions have 
odd weights. This by no means implies that the great majority of totally- 
irreducible functions have odd weights for growing K. Indeed from (34) 



So, half of all the i^-Boolean functions are totally-irreducible with an odd 
weight. Now, for the complement %k (K) \Zk (the totally-irreducible func- 
tions with even weights): 

Since (3k given by (15b) has the asymptotic behavior (16) for K >> 1; 



So, with respect to the normalized counting measure, almost half of the 
totally-irreducible ^-Boolean function have an even weight. 

Normalizing according to 2 2 * -1 , totally- irreducible ^-Boolean functions 
with odd weights distribute by a binomial distribution with mean (tju) = 
2 2 ~ 1 [ as needed for consistency with (33)], and a standard deviation E 2 = 
2^~ 2 ; so the relative error E/ (to) = 2~ K / 2 , decays exponentially to zero with 
growing K. For K ^> 1 the binomial becomes asymptotically a gaussian 
distribution, with the very same moments. 



#3^ = 2 
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5. Conclusions 

In Ref. [9] a classification of Boolean functions, in terms of their irreducible- 
degree of connectivity, was constructed allowing to calculate exact, and 
asymptotic behavior of iVi^-Kauffman networks. In this work we have devel- 
oped further the mathematical structure and consequences of this classifica- 
tion with the aid of the well known ring-isomorphism (17) i.e. 2 n = VVt 14 . 
In Theorem /we found an operational formula (26) for the irreducible-degree 
A (bx) of a fT-Boolean function that allows concrete analytical manipulations 
of it for calculational purposes. 

In Theorem II, we established a ring-isomorphism between the set 
D\ K . . . , ix} of reducible fT-Boolean functions on the indexes . . . , i\} C 
[K] and the double power set V 2 ([K] \ {ii, . . . , This fundamental theo- 
rem allows to use the rich ring-isomorphic structure for analytical calculations 
over Boolean functions. In particular, we have used it to prove Theorem III 
which gives a formula for qk (A, oj): the number of i^-Boolean functions with 
irreducible-degree A and weight oj. The value of qk (A, oj) turns out to be 
of fundamental importance in calculations of quantities which give informa- 
tion about the dynamics of iVif-Kauffman networks, as: i) The probability 
of extracting a i^-Boolean function with irreducible-degree A, through the 
distribution II (u) Eq. (8). This quantity is an ingredient in the study of the 
stability of NK -K&u&m&n networks as see in Appendix A Eq. (A2). ii) The 
study of the phase transition of the A^-Kauffman networks done in Ref. [10] 
(see Appendix B for a summary). 

We are using this framework further as a tool for the calculation of the 
dynamical properties of A^-Kauffman networks with promising results. 
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Appendix A: About the Canalizing Boolean Functions 

For clarification purposes, we state now in the mathematical framework of 
this article the concept of a canalizing i^-Boolean function, which was first 
introduced by Kauffman 3 ' 12 . 
Definition: 

A bx '■ %2 — > Z 2 Boolean function is said to be canalizing iff 3 i G [K] , 
and £, r G Z 2 such that 

Si = £ =4> b K (Si, . . . , Si-i, £, S i+ i, . . . , S K ) — r, 

V Sj G Z 2 with j G [K] \ {i}. 

As examples for K = 2: in Table 1 all Boolean functions, but 7 and 10, 
are canalizing. In more detail: function 15 = V also known as the OR 
function, is canalizing since if one their two arguments (Si or S 2 ) is set to 
"1", then OR = 1 whichever the value of the other argument. We see that 
this definition divides E K in two disjoint subsets (the canalizing, and the 
no-canalizing sets). 

N.B. Canalization and irreducibility are different properties of Boolean 
functions, neither of them containing the other. While canalizing functions 
are partially responsible of the frozen dynamics of Ai^-Kauffman networks 
since they generate forcing structures that tend to trap the dynamics into 
short loops 3 ; they do not offer a full understanding of their dynamics. This 
requires a different and shifter information of the behavior of bx as a func- 
tion of their arguments; as irreducibility does 9 . Irreducibility appears in a 
natural way in the calculation of important quantities that govern the dy- 
namics of Aif -Kauffman networks 9 ~ n > 16 . We quote three examples from our 
researches: 

i) The calculation of the mean number of Ai^-Kauffman networks mapped 
to the same functional graph $(N,K), was done in the Appendix of 
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Ref. [9] equations (R18), and (R19); with the result 
&(N,K) = {l-<p(N,K)y N , 



(R18) 



where 



<p(N,K) 



Eto^(A)[(^)-l] 



(#19) 



As we see, this important quantity, is expressed as a series generated by 
Pk (A) (the number of Boolean functions with irreducible-degree A). So, 
decomposition (11) with their cardinalities (15) contain the information 
of the dynamical diversity of A^-Kauffman networks. Furthermore, 
each term of the series becomes smaller with growing A; which yields, 
for N 1, and K ~ O (In In N), the asymptotic expression 9 ' n 



Using this result into (R18) and looking for a value K c (of K) such that 
•&(N,K C ) = 1/2, it is obtained that a critical connectivity K c exists, 
such that: 

a) It is given by 



b) For K < K c {) (N, K) >> 1, so many A^-Kauffman networks are 
mapped into the same functional graph. 

c) For K > K c ; d (N, K) I. So in this case, almost any NK- 
Kauffman network is mapped to a different functional graph. 

See Refs. [9,11] for details. 

ii) The calculation for the probability P (A); that an A^-Kauffman net- 
work [Appendix B Eq. (Bl)] remains invariant against a change in 
one of their ^-connection functions C*^ was reported in Eqs. (23) of 
Ref. [9] (as one of the main ingredients for final result) obtaining, 




K c log 2 log 2 




K 



K\ (N-X)\ 
N\(K-\y 



P [b K e % K (A)] , 



(Al) 



A=0 
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where P [b K E %k (A)] is the probability to extract a function with 
irreducible-degree A. Once again we see other important quantity that 
gauges the dynamical behavior of iViT-Kauffman networks and is ex- 
pressed by a series of quantities that are function of decomposition (11) 
of E K in their irreducible subsets %k (A). 

Note that P [bx G (A)] is analytically expressed by the use of (8) 
and (31) giving 

2 K 

P [b K e % K (A)] = sk (A, u) n (u) . (A2) 

ui=0 

In Ref. [9], however, the problem was focused in the asymptotic behav- 
ior of (Al) for N > 1, and K ~ O (1), thus obtaining 

P(A)*P [b K G I K (0)] + O (£) = P 2K + (1 - p) 2 " + O (1) 

without requiring the full calculation of (A2). 

Hi) The transition curve for the dynamics of Ai^-Kauffman networks in a 
mean field treatment is shown to depend in the average connectivity, 
which is a function of A {bx) and not of canalization. See Appendix B, 
where the transition curve (B2) is obtained as the average of A (bx) 
weighted by (8) and (31); see also Ref. [10] for a more detailed study 
of the mean field treatment. 
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Appendix B: Mean Field Dynamics of A^-Kauffman Networks 



We summarize the mean field approach, corrected for irreducibility, to 
study the dynamics of the iVi^-Kauffman networks. A detailed study should 
be consulted in Ref. [10]. 

iVif-Kauffman networks, are Boolean endomorphisms / : Z^ — > Z 2 of 
the form 

Z% ^ Zf ^ Z 2 i = l,...,N, 

where the connection function C*® : Z 2 — » Zf- cuts N — K oi the N 
Boolean variables Sj, j = 1, . . . , N, so C*® (Si, . . . , Sn) = (S^, . . . , S i[C ), 
with {ii, . . . ,i K } C [N] (i a G [N], a G [K]) being whichever of the (^) sub- 
sets with cardinality K of [N\. The dynamic is defined by the synchronous 
iterations 

S- i (* + l) = 6goC^ (i) (S(*)), i = l,...,iV. (51) 

The 6k Boolean functions are extracted randomly according to the probabil- 
ity distribution (8), while the connection functions C*^ are extracted with 
equiprobability from the (^) possible ones. 

The problem, is to study how is the way in which the dynamics of (Bl) 
behaves as a function of their defining parameters, which are N, K, and the 
bias p of the probability distribution (8) with which the bx functions are 
extracted. 

A way to observe how the dynamics of (Bl) evolves is to see the behavior 
of the Hamming distance of two nearby states S, S' G Z^ 

TV 



d H (s, s') = + si\ 



i=i 



for asymptotically big values of N. The dynamic generated by (Bl) is deter- 

(i) *(i) 

ministic, but the construction of the functions b K and C K which determines 
the endomorphism is done randomly. This allows to do a statistical treat- 
ment of the dynamics and make a mean field approximation 2 ' 4 > 10 . For that: 
let us see, that due to the randomness of the construction of (Bl), each site 
% at t — 0, such that Si(0) ^ 5^(0), will affect on average K sites; each one of 
them is going to be, also, the argument of a bx function at the next iteration 
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t — 1. So, at the next step, site i will contribute on average to the Hamming 
distance by the factor 

*i = (P e (b K ) X(b K )) i , 

with the average taken with respect to the bf) that contribute to the i- 
site, and P c (bK) being the probability that each bf) changes its output, due 
that one of their arguments has changed, (which is explicitly calculated in 
Ref. [10]). If the system starts at an initial Hamming distance d#(0) = 
dtf (S(0),S'(0)) such that 1 < d H (0) < N, for N > 1 we may apply the 
central limit theorem and take the average over the iV-sites 

<*> = -^ ^i = (Pc(b K ) X(b K )). 
i=i 

Now, at t — 1 Hamming distance will grow (or decay) on average by 

dH(l) wdn(0) (^c(6at) A(&*)> , 
or more generally, while the condition 1 <C d H (t) <C N is fulfilled 

d H (t + l)nd H (t) (P c (b K ) X(b K )). 
Solving this difference equation we obtain 

d H {t) ~ d H (0) (Pc(b K ) X(b K )y. 

So depending on whether (P c {Pk) X(bx)) is greater or lower than 1, Hamming 
distance will grow or decay exponentially; with the equation 

A(K,p) = (P c (b K ) X(b K )) = l 

signaling the phase transition frontier. A (K, p) may be calculated explicitly 
through the use of (31) obtaining the equation 

2 K K 

A (K,p) = ^2 n (w) P c (w) ]T A q k (A, w) (52) 

w=0 A=0 

which is explicitly calculated in Ref. [10], obtaining; 
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A(K,p) = 2Kp(l-p)x 



x {l - 2p (1-p) [1 - 2p (1 -p)] 2 ^ 1 - 2 } = 1, (53) 

for the phase transition curve. Equation (B3) is an improvement of a result 
previously obtained in 1986 by Derrida & Stauffer Ref. [15], for the mean 
field treatment, where reducibility of Boolean functions was not taken into 
account; obtaining the result 

A(K,p)=2Kp(l-p) = l 

for the transition curve. See Ref. [10] for more details. 
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